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ABSTRACT 

We have determined the strong coupling as from a comprehensive study of 
energy-energy correlations {EEC) and their asymmetry (AEEC) in hadronic 
decays of bosons collected by the SLD experiment at SLAC. The data were 
compared with all four available predictions of QCD calculated up to O^a"^) in 
perturbation theory, and also with a resummed calculation matched to all four 
of these calculations. We find large discrepancies between as values extracted 
from the different 0{a'^) calculations. We also find a large renormalization scale 
ambiguity in as determined from the EEC using the 0{a'^) calculations; this 
ambiguity is reduced in the case of the AEEC, and is very small when the matched 
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calculations are used. Averaging over all calculations, and over the EEC and 
AEEC results, we obtain as(M|) = 0.1241 Jgg|(ea;p.) ± O.Q09{theory). 
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1. I. INTRODUCTION 

Measurement of the strong coupling as in various hard processes and at 
different hard scales is one of the crucial tests of Quantum Chromodynamics 
(QCD) [1]. In e"'"e~ annihilation ag may be determined from inclusive measures of 
the topology of hadronic events. We have previously determined as from the rate 
of multijet events in hadronic decays of bosons collected by the SLD experiment 
at SLAG [2]. We found as = 0.118 ± 0.002(stot) ± 0.003(s?/st.) ± 0.010(i/ieory), 
where the dominant uncertainty arises from uncalculated higher order contributions 
in perturbation theory. Here we present complementary measurements of as 
using energy-energy correlations {EEC) and the asymmetry of energy-energy 
correlations {AEEC) [3]. These are inclusive two-particle correlations that can 
be used to probe the structure of hadronic events in more detail than the event 
topology variables and can be calculated perturbatively in QCD. Comparison of as 
determined in this way with that measured from event topology variables provides 
a significant consistency check of the validity of perturbative QCD. 

The EEC is defined as the normalized energy-weighted sum over all pairs of 
particles whose opening angles Xij lie between x ~ and x + 

event ^ J .7^, E^^^ 

A 2 

where x is an opening angle to be studied for the correlations; Ax is a bin width; 
Ei and Ej are the energies of particles i and j and E^is is the sum of the energies 
of all particles in the event. In the central region, x ^ 90°) the shape of the EEC 
is determined by hard gluon emission; hadronization contributions are expected to 
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be large in the collinear and back-to-back regions, x ~ 0° and 180° respectively. 
The asymmetry of the EEC is defined as 



AEECix) = EECin - x) - EEC{x) ■ (2) 



Perturbative QCD calculations of the EEC were first performed up to 0{ag) 
in 1983 by Richards, Stirhng and EUis (RSE) [4]. Since then similar calculations 
have been performed by Ali and Barreiro (AB) [5], Falck and Kramer (FK) [6], 
and Kunszt and Nason (KN) [7]. These calculations, valid in the central region, 
have the general form 

EEC{x) = ^^MX) + (^^) ' [Mx)27ibo Hi^^s) + B{x)] , (3) 

where, to the same order in perturbation theory, as{fJ>^) is related to the QCD 
scale A.j^ by [8] 



MS' 



1 



5i ln[ln(/.VA|^)] 
bl ln(/i2/A.2 ^ 



MS' 



(4) 



ji is the renormalization scale, often expressed in terms of the factor / = yU^/s; 
^/s is the center-of-mass energy of the experiment; 6o = (33 — 2nj)/127r; 
bl — (153 — 19nj^)/297r2; and is the number of active quark fiavors. Here 
we have assumed the definition of A;^ for five active fiavors. The first order 
coefficients A{x) can be calculated analytically and the second order coefficients 
B{x) are calculated numerically. The main difference between the four theoretical 
calculations mentioned above is in the method used to treat the singularities 
appearing in the second order coefficients in the angular regions x ~ 0° and ~ 180°. 



It is important to note that 0{a1) perturbative QCD calculations do not 
specify the /i value to be used for any physical observable, although this scale 
ambiguity will presumably vanish if the calculation is done to all orders in 
perturbation theory. Large scale ambiguities in determinations of as using such 
calculations for the EEC and event variables have been reported [2,9-11]. There 
is an indication that the AEEC may be less sensitive to higher-order perturbative 
QCD contributions than the EEC [4], and may therefore be expected to be less 
sensitive to changes of renormalization scale. One also expects a priori that 
non-perturbative hadronization effects will tend to cancel in the AEEC. 

Furthermore, discrepancies between the four 0{al) calculations of the EEC 
have been reported [11]. These discrepancies between calculations each supposedly 
complete to 0{a1) are not understood, but may be interpreted as an indication 
that not all 0{a^g) terms have been included in some or all of the calculations. In 
the absence of further information we assume that all calculations are equally valid 
and use them all in our determination of as, taking the spread in ag values as an 
indication of calculation uncertainty. 

Recently progress has been made in perturbative QCD in the form of 
'resummed' calculations of certain event shape measures in e^e~ annihilation [12]. 
These techniques have been used [13,14] to calculate the EEC in the back-to-back 
region, where a large contribution from soft and coUinear gluon radiation appears, 
by exponentiating the leading and next-to-leading order terms in L = ln(l/y), 
where y = ^^^^^^ , up to all orders^ in ag ■ Within this formalism the cumulative 
EEC can be written [13] 

f Earlier work on a perturbative evaluation of the leading and next-to-leading order terms in 
L of the EEC was performed up to two-loop level by Collins and Soper [15]. 
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^ oo oo 

REEciy) = / EEC{y')dy' = (l + a^C„) a,) + (5) 

Q n=l n=l 

where 

oo n+1 

ln(E(y,a,)) = 5] Gn^ln-(1/|/) 

n=l m=l 
= ( ^121/^ + Gi\L ) as 

+ ( G2ZL^ + 6-22^^ + ^21^ ) Oil 
+ ••• 

^ Lgi{asL) ^ + ^ g2{asL) ^ + ^ aggsjasL) H j (6) 

Leading Next-to-leading Subleading 

Cn and Gnm are constants; Yn{y) are functions which vanish as y — > 0; and the 
functions gi and g2 are the sums of the leading and next-to-leading logarithms in 
L, respectively. Except for the term proportional to G21, which can be estimated 
from the 0{ag) calculations, the subleading terms have not been calculated. An 
analytic evaluation of the EEC singularity structure [13] gives an approximate 
simplified form: 

Rmny) = (1 + Cias) exp{Lgi{asL) + g2{asL)). (7) 

This resummed calculation can be combined with each of the 0{al) calculations, 
also expressed in the cumulative form: 

= 1 + A{y)as + B{y)al (8) 
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where A and B are the cumulative forms of A and B in Eq. (3). Subtracting 
double-counted terms yields the matched form: 

ReecHv) = REEcPiy) + R'^EEciy) - (Ci + GuL + Gr2L')as (9) 
- [G22L'' + G23L^ + {GuL + G12L') (Ci + ^ (GiiL + daL^))] 0.2. 

This procedure is called the ' R-matching^ scheme.^ Equation (9) can be 
differentiated with respect to x ^^'^ compared with the data. A 'modified 
R-matching' scheme is also proposed in Ref. [13], whereby the subleading term 
proportional to G21 at 0{a^g) (Eq. (6)) is included in the argument of the 
exponential in Eq. (9). It can be argued that modified R-matching is preferred 
theoretically [17] because the extra subleading term at 0{al) is explicitly 
exponentiated. 

In this analysis we attempt to be as comprehensive as possible and compare 
all four 0{al) calculations with data in the central region to determine as- We 
then combine the resummed calculation with all four of the fixed order calculations, 
using both the R-matching and modified R-matching schemes, to make an improved 
determination of as- We study in particular the dependence of as on the QCD 
renormalization scale, which is expected a priori to be weaker when the resummed 
calculation is used. We compare our results with as measured from our jet rates 
analysis [2] and with measurements from LEP [9-11,16,18-22]. 



f Another scheme, HnR-Tnatching\ has also been proposed [12] to combine the resummed and 
0{al) calculations for event shape variables, but it cannot be applied reliably to the EEC 
[16]. 



7 



2. II. THE SLD AND EVENT SELECTION 

The data used in this analysis were recorded in 1992 by the SLC Large Detector 
(SLD) [23] from electron-positron annihilation events at the resonance produced 
by the SLAC Linear Collider (SLC). A subset of the components of the SLD is used 
for this analysis. Charged particles are tracked in the vertex detector (VXD) [24] 
and in the Central Drift Chamber (CDC) [23]. Momentum measurement is 
provided by a uniform axial magnetic field of 0.6 T. Particle energies are measured 
in the Liquid Argon Calorimeter (LAC) [25], which contains both electromagnetic 
and hadronic sections, and in the Warm Iron Calorimeter [26]. In this analysis the 
calorimeters are used only for triggering and event selection. 

Three triggers were used for hadronic events, one requiring a total LAC 
electromagnetic energy greater than 30 GeV, another requiring at least two 
well-separated tracks in the CDC, and a third requiring at least 8 GeV in the 
LAC and one track in the CDC. A selection of hadronic events was then made by 
two independent methods, one based on the topology of energy depositions in the 
calorimeters, the other on the number and topology of charged tracks measured in 
the CDC. 

The analysis presented here used charged tracks measured in the CDC and 
VXD. A set of cuts was apphed to the data to select well- measured tracks and 
events well-contained within the detector acceptance. Tracks were required to 
have: (1) a distance from the measured interaction point, at the point of closest 
approach, of less than 5 cm in the direction transverse to the beam axis and 10 cm 
along the beam axis; (2) a polar angle Q with respect to the beam axis within 
jcos^l < 0.80; (3) a momentum transverse to the beam axis of pj^ > 150 MeV/c. 
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Events were required to have: (1) a minimum of five such tracks; (2) a total visible 
energy calculated from the selected tracks of at least 20 GeV; (3) a thrust axis [27] 
direction within |cos^^r| < 0.71. Prom our 1992 data sample 6476 events survived 
these cuts. The efficiency for selecting hadronic events satisfying the |cos6't| cut 
is estimated to be above 96%. The background in the selected event sample is 
estimated to be 0.3±0.1%, dominated by Z'^ — > t~^t~ events. Distributions of 
single particle and event topology observables in the selected events are found to 
be well described by Monte Carlo models of hadronic boson decays [28,29] 
combined with a simulation [30] of the SLD. 

3. III. MEASUREMENT OF EEC AND AEEC 

The EEC was calculated using all pairs of selected charged tracks, assigning 
each the charged pion mass. Figure 1 shows the measured EEC and AEEC. The 
bin width was chosen to be 3.6°, which is much larger than the two-particle angular 
resolution of the detector, so as to minimize bin-to-bin migration effects in the data 
correction procedure to be described later. Also shown in Fig. 1 are comparisons 
with the JETSET [28] and HERWIG [29] Monte Carlo programs which simulate 
the hadronic decays of Z° bosons, combined with a simulation of the SLD, and 
the same event reconstruction and selection procedures as applied to the data. For 
each program 100,000 events were generated. For JETSET we used parameter 
values determined by the TASSO Collaboration at ^/s — 35 GeV [31], which have 
been found to be in good agreement with Z° data [32] . For HERWIG we used the 
default parameter values which were derived from comparisons with LEP data [33] . 
Both simulations reproduce the general features of the data. 
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Wc used both JETSET and HERWIG simulations to correct our data for 
the effects of initial state photon radiation, detector acceptance and resolution, 
interactions and decays in the detector, analysis cuts, and unmeasured neutral 
particles. The measured EEC, EECmeasiXi)^ was corrected to the hadron level, 
EEChadroniXi), by applying bin-by-bin correction factors, C£f^(xj): 

EEChadroniXi) = Cgf^ixi) X EECmeasiXi), (10) 

where the correction factors were calculated by comparing Monte Carlo results 
before and after detector simulation: 

^^'^detectoAX^! 

where EEC^^^^^^Xi) represents the histogram content at bin Xi of the EEC 
obtained from the charged particles of the reconstructed Monte Carlo events, and 
EECf^^^g^{Xi) represents that from all the stable particles with hfetimes greater 
than 3 x 10~^° s in the generated Monte Carlo events with no initial state photon 
radiation. The corrected AEEC was then derived from the corrected EEC. 
Systematic errors in this correction procedure will be discussed later. 

We corrected the data further for the effects of hadronization using both 
JETSET and HERWIG. The parton-level EEC is defined by 

EECpartoniXi) = C'fm^(Xi) X EEChadroniXi), (12) 

where the correction factor C^^^{xi) is defined by 

riEEC(^, \ _ ^^^partoniXi) , , 

'-/raff \Xi) - j^j^^MC (^.y ^^-^f 
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where EECffJ;l^^{xi) is the hadron- level EEC as defined in Eq. (11), and 
EEC^^f.^{Xi) is the content of bin Xi constructed at the parton level? The 
corrected AEEC was again derived from the corrected EEC. 

The parton-level EEC and AEEC are shown in Figs. 2(a) and 3(a) 
respectively, where JETSET was used for both the detector and hadronization 
corrections. The correction factors for the EEC for detector effects and 
hadronization effects are shown in Figs. 2(b) and (c) respectively. Figures 3(b) 
and (c) show the corresponding effective correction factors for the AEEC, which 

are defined as Cj^^^ [xi] = AEEC^^g^^{xi)/AEEC^^{xi) and Cf^^^f{xi) = 
AEEC^^^t^{Xi)/AEEC^^g^{Xi) respectively. Also shown in Figs. 2(b) and (c), 
and Figs. 3(b) and (c) are the correction factors calculated using HERWIG. It is 
interesting to note that JETSET and HERWIG give slightly different correction 
factors in the central region of the EEC (Figs. 2(b) and (c)), although both 
describe the detector-level data well (Fig. 1(a)). We have found that this is due to 
differences in the relative production of charged and neutral particles, which also 
lead to the small differences seen in Fig. 1. However, the overall corrections to 
the parton level are found to be very similar in the central region. The differences 
between the JETSET and HERWIG correction factors will be discussed in the 
section on systematic errors (Section IV). It is also interesting that, despite the 
a priori expectation that hadronization effects should cancel in the AEEC, the 
magnitude of the deviation of the hadronization correction factors from unity is 
comparable for the AEEC and EEC. 

t The parton level is defined in JETSET by the parameter Qq, which was set to 1 GeV, and 
in HERWIG by the parameter Qg, which was set to 0.75 GeV. 
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The statistical errors on the EEC and AEEC have strong bin-to-bin 
correlations because each track in each event contributes to several bins. The 
statistical error in each bin was estimated by taking the rms of the contents of 
that bin over ten Monte Carlo samples, each with the same number of events as 
the data sample. These errors were used in the fits below. 

4. IV. DETERMINATION OF ag 

A. Determination of as using 0{ag) calculations 

We first determine ag by comparing the four 0{ag) QCD calculations with 
the parton-level EEC and AEEC derived from the data. Each calculation was 
fitted to the fully corrected measured EEC and AEEC by minimizing under 
variation of A-^^ for fixed renormalization scale factor /. The fits were restricted 
to the angular region 36° < x < 154.8° for the EEC and 21.6° <x< 79.2° for 
the AEEC, the ranges being chosen in such a way that (1) the uncertainty on the 
hadronization corrections is smaller than 30%, (see Figs. 2(c) and 3(c)), (2) the 
perturbation series for the EEC and the AEEC remain reasonably convergent 
for all four QCD calculations [34], and (3) the fit quality at / = 1 is reasonable, 
Xdof ^ 5, where x'dof P^^^ degree of freedom of the fit. 

Figures 2(a) and 3(a) show the results of the four 0{ag) QCD fits with 
renormahzation scale factor / = 1, to the EEC and AEEC respectively. 
Figures 4(a) and 5(a) show as{M'^). derived from the fitted A-^^ for diff'erent 
values of /, using the EEC and AEEC respectively. The corresponding fit quality 
Xdof is shown in Figs. 4(b) and 5(b). In order to estimate the statistical errors 
on as we made use of the previously generated ten sets of Monte Carlo events. 
We performed the same fitting procedure to the EEC and the AEEC for each of 
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these sets and took the rms deviation of the ten as values thus determined to be 
the statistical error of the fitted ctg (see Figs. 4(a) and 5(a)). 

Prom Figs. 4 and 5 several features of the fit results are common to each QCD 
calculation: (1) q;s(M|) depends strongly on /, the dependence being stronger for 
the EEC than for the AEEC; (2) at low / the fit quahty deteriorates rapidly, and 
neither nor its error can be interpreted meaningfully; (3) the fits are relatively 
good in the scale range / ^ 2 x IQ-^ for the EEC and / ^ 5 x IQ-^ for the AEEC, 
and there is no strong preference for a particular scale. Similar features were seen 
in our ctg measurement from jet rates [2] . It should be noted that a comparison of 
as results from the various calculations at any particular scale / reveals systematic 
differences of up to 20% in magnitude,^ although for the AEEC the three more 
recent calculations are in reasonable agreement for / > 0.01. 

Figures 4 and 5 represent a complete summary of our results. However, in order 
to quote a single value of as from each QCD calculation we adopt the following 
procedure. We consider the ranges 0.002 < / < 4 for the EEC and 0.01 < / < 4 
for the AEEC. The lower bounds ensure that the perturbation series for the EEC 
and AEEC remain reasonably convergent for all four QCD calculations [34] and 
that the fit quality is reasonable. The upper bound restricts to a reasonable 
physical region, < 2^/s. We take the extrema of the as values in these ranges 
to define the scale uncertainty. Table 1 summarizes the measured ttg, defined as 
the midpoint of as between the extrema, and the scale uncertainty, defined as the 
difference between the extrema and the midpoint value. 



f This has also been observed previously [11]. 
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The experimental systematic errors, which arise from uncertainties in modehng 
the acceptance, efficiency, and resolution of the detector, are summarized in 
Table 2. The largest contribution is from the understanding of the effects of our 
track and event selection cuts, which we evaluate by varying the cuts over wide 
ranges. We also varied the tracking efficiency and resolution by large amounts in 
our Monte Carlo simulations to account for any possible biases. In addition, the 
change in ag is negligible when bins at either end of the fit range are removed. 
Effects due to limited Monte Carlo statistics are relatively small compared with 
the other errors. Total experimental systematic errors on as are estimated to be 
+2:SS3 for the EEC and ±0.004 for the AEEC. 

The theoretical uncertainties in parton production and hadronization were 
estimated by recalculating the hadronization correction factors using JETSET with 
values of Qq in the range 0.5-4.0 GeV,^ and by using HERWIG, which contains 
a different hadronization model. The uncertainties in the correction factors are 
shown in Figs. 2(c) {EEC) and 3(c) [AEEC). The estimated hadronization 
uncertainties in as are ±0.002 for the EEC and AEEC^ 

The theoretical uncertainties from the numerical precision of the coefficients 
A{x) and B{x) in Eq. (3) are found to be negligible for all four QCD calculations. 

Averaging over all four 0{al) QCD calculations to quote a single as value from 
each of the EEC and AEEC analyses, we calculated the mean and rms deviation 
of the four as values at each / in the ranges shown in Figs. 4 and 5, respectively. 
The central value of ag was taken as the midpoint of ag in that / range; the 

f Differences observed for changes of the other main fragmentation parameters in JETSET, 

Uq, a and 6, are small and negligible compared to those of Qo- 
I These should be compared with overall shifts in q;s(M|) of +0.014 {EEC) and —0.015 

{AEEC) due to the hadronization corrections. 
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calculation uncertainty was taken as the rms at the central value; and the scale 
uncertainty was taken as the difference between the central value and the extrema. 
We obtain: 

EEC : as{M^) = 0.127 ±Q.QQl{stat.)toS{syst.) ± QM3{theory) 

AEEC : aa{Ml) = 0.116 ± 0.002(stot.) ± 0.004(s|/st.) ± 0.006(t/ieor|/). 

The total theoretical uncertainty is the sum in quadrature of contributions from 
hadronization (±0.002 for EEC, +J:oo2 AEEC), calculation (± 0.006 for 
EEC, ± 0.004 for AEEC), and scale (± 0.011 for EEC, ± 0.003 for AEEC) 
uncertainties. These as values are in agreement with our measurement from 
jet rates [2]. They arc also in agreement within experimental errors with as 
measurements from the EEC and AEEC by other groups [9-11,18]. A large 
theoretical uncertainty, dominated by the scale uncertainty, was also found in our 
as measurement from jet rates [2] and is related to uncalculated higher order terms 
in QCD perturbation theory. The scale uncertainty for the AEEC is smaller than 
that for the EEC, perhaps indicating that the contribution from uncalculated 
higher order terms is smaller in the former case as discussed in Section I. Finally, 
the discrepancy in as between the four calculations is a significant contribution to 
the total uncertainty for the EEC, and the dominant contribution for the AEEC. 

B. Comparison of 0(0;^) and (9(Q;g)+resummed calculations 

We now compare the 0{a^g) calculations with the matched 0(Q;g)+resummed 
calculations. For illustration we use the KN 0{a^g) calculation and the modified 
R-matching scheme. For this purpose we perform fits to the EEC data within 
the angular range 90° < x ^ 154.8°, where the lower limit is the kinematic limit 
for the resummed calculation and the upper limit is the same as in the previous 
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section. Figure 6 shows the results of 0(0^) and matched fits, with / = 1, to 
the parton-level corrected EEC. At / = 1, as{Ml) = 0.126 ± 0.002(stoi.) with 
;^2^^ ^ 1.86 for the 0{al) fit, and as(M|) = 0.133 ±0.002(stat.) with xlof = O-'^l 
for the matched fit. However, the 0{a1) calculation cannot describe the large 
angle region x ^ 150°, while the matched calculation describes the data even up 
to X ~ 170°. We show in addition in Fig. 6 a fit of the 0{a1) calculation with 
/ = 0.0033, which results in q;s(M|) = 0.114 ±0.001(stat) with xlof = With 
this small scale choice the 0{a1) calculation describes the data up to x ~ 165°. 
This illustrates that an 0{al) calculation used with a small scale is able to 
reproduce, to some extent, the effects of the resummed logarithms in the matched 
calculation. However, the small scale choice results in a smaller fitted value of 
q;s(M|) than the choice f — 1. 

Figure 7 shows the resulting q;s(M|) and Xdof values for fits of the and 
matched calculations at various values of /. It is apparent that good fits to the 
data can only be obtained for / ^ 10~^ (0(0^)) or / ^ 0.2 (matched).^ Within 
the range / ^ 0.2 the scale dependence of the matched calculation is significantly 
smaller than that of the 0{a'^) calculation, confirming the a priori expectation 
discussed in Section I. 

C. Determination of ag using 0(Q;g)+resummed calculations 

In this section we apply both the R-matching and modified R-matching schemes 
(see Section I) to combine the resummed calculation with each of the four 0{al) 
calculations of the EEC, and perform fits to the data. Note that such matched 

f Similar results were obtained in our measurement of as from jet rates [2]. 
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calculations cannot be applied to the AEEC because no resummed calculation is 
available in the forward region, 0° < x < 90°. 

Figure 8 shows the results of the KN C(ci;g)+resummed fits to the EEC, 
applying R-matching and modified R-matching with renormalization scale factor 
f — 1- The fit range was restricted to the angular region 90° < x < 154.8° for 
R-matching and 90° < x < 162° for modified R-matching according to the criteria 
discussed in Section IV.A. The R-matching case does not describe the data in the 
large angle region, x ^ 155°. Figure 9 shows the resulting as{M^) and Xdof values 
as a function of / for modified R-matching with all four 0{al) calculations. The 
most significant features are that / ^ 0.2 is needed to fit the data, and ctg varies 
slowly with /, resulting in a significantly smaller scale uncertainty than in the 
0{a^) case, as expected. Once again, large differences in ag are apparent between 
the matched RSE, AB, FK and KN calculations. The results for R-matching, 
shown in Fig. 10, are qualitatively similar, but are found to have a slightly larger 
variation of ag with / than that from modified R-matching, probably due to the 
neglect of the sub-leading term in the exponentiation (Section I). 

To quote a single value of ag for each matched calculation we adopt the same 
procedure as in the fits but use the range 1/4 < / < 4, where the lower 

17 



bound ensures that fit qualities remain reasonable. This is the same / range 
that we considered in our ctg measurement from jet rates using resummed 
calculations [2]. The results are shown in Table 1. We note that in each case the 
matched fit yields systematically larger as values than the 0{ag) fit, and that the 
modified R-matching results are systematically larger than the R-matching results. 
Taking the average over all four calculations as before, and averaging over both 
matching schemes, we obtain 

as{Ml) = 0.130 ±0.002{stat.)^^^^^^{syst.) ± 0.007 (theory), 

where the total theoretical uncertainty is the sum in quadrature of contributions 
from hadronization (± 0.002), calculation (± 0.005), scale (± 0.002), and 
matching (± 0.004) uncertainties. Here we included a matching uncertainty from 
the difference in as values between R-matching and modified R-matching as a 
contribution to the total theoretical uncertainty. It is interesting to note that 
the dominant uncertainty in this case arises from the discrepancies between the 
four 0{ag) calculations, the scale uncertainty being even smaller than the total 
experimental error. Resolution of this discrepancy between calculations would 
potentially yield, therefore, a precise measurement of as in e+e" annihilation at 
the resonance. The above result is in agreement with that from the 0{a'^) fits 

f For R-matching, scale factors as small as / ^ 0.06 arc allowed by the quality of the fits 
to the data (Fig. 10(b)). However, Fig. 10(a) shows that the change in a,,(M|) between 
/ = 1/4 and / — 0.06 is at most 0.002 in magnitude (FK calculation). The change in 
the mean as{M^), averaging over all four calculations, is smaller than 0.0004, which is 
negligible compared with statistical errors. Therefore, extending the / range down to 0.06 
does not affect our final result, and we use the range 1/4 < / < 4 for consistency with 
our jet rates analysis [2] . 
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to the EEC within experimental errors, and with the result from the 0{ag) fits to 
the AEEC within the combined experimental errors and theoretical uncertainties. 

5. V. SUMMARY 

We have measured energy-energy correlations and their asymmetry in hadronic 
decays and compared our corrected data with four 0{a1) perturbative QCD 
calculations in order to extract as- We have also combined these calculations with 
a resummed calculation using two matching schemes and extracted ag from the 
EEC. We obtained 

as{Ml) = Q.m^^^^^\exp.) ± Qm2>{theory) from EEC {0{al)), 

as{Ml) = 0.116 ± 0.005(ea;p.) ± 0.006(t/ieor|/) from AEEC (C(a^)), 

as{Ml) = 0.130 ^Jj'|J°^(exp.) ± 0.007(i/ieory) from EEC (C>(q;2) + resummed), 

where the experimental error is the sum in quadrature of the statistical and 
systematic errors. Note that all experimental errors and theoretical uncertainties 
are highly correlated. These as values correspond to A-^^^^ ~ 360MeV, 200MeV 
and 420MeV, respectively. The renormahzation scale ambiguity dominates the 
uncertainty in the first case, but in the second and third cases it is dominated 
by the discrepancy between the four 0{al) calculations. In the third case 
there is also a large uncertainty from matching the 0{a1) and resummed 
calculations. The renormahzation scale dependence of the C(tts) results, as 
well as the reduction of renormahzation scale sensitivity and the increase in 
measured ag with the application of resummed calculations, is similar to that 
reported in our measurement of ag from jet rates [2]. The results using 0{a1) 
and resummed-|-C(Q;g) calculations are also consistent within experimental errors 
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with similar analyses by the LEP experiments [9-11,16,18-22]. Taking an 
unweighted average over all three results we obtain: q;s(M|) = 0.124+jJ;jjg(ea;p.) ± 
0.009{theory), corresponding to Aj^ = 317ttl{exp.) ± UA{theory) MeV. 
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Table 1: as{M'^) and scale uncertainties from 0{a^g) QCD fits to the EEC and 
the AEEC, and from C»(a2)+resummed fits to the EEC. 



QCD 
calculation 


El 

as{Ml) 


0{al 

Scale 

uncertainty 


) fits 

AE 

as{Ml) 


EC 

Scale 

uncertainty 


( 

mod. R- 

a,(M|) 


D{a1)-\-Yes\ 
EE 
matching 

Scale 

uncertainty 


immed fil 

;c 

R-ma 

as{Ml) 


s 

tching 

Scale 

uncertainty 


AB 
FK 
KN 

RSE 


0.132 
0.119 
0.125 
0.133 


± 0.011 
± 0.013 
± 0.012 
± 0.011 


0.114 
0.113 
0.114 

0.124 


± 0.004 
± 0.003 
± 0.004 
± 0.008 


0.138 
0.127 
0.133 
0.140 


± 0.001 
± 0.002 
± 0.001 
± 0.001 


0.130 
0.120 
0.124 
0.130 


± 0.003 
± 0.004 
± 0.004 
± 0.003 



Table 2: Summary of estimated experimental systematic errors on ag (M|) for 
the EEC and AEEC. 



Source 


EEC 


AEEC 


Tracking 






Track selection 


+0.0006 

-0.0011 


+0.0006 

-0.0011 


Tracking efficiency 


+0.0006 


+0.0013 


Momentum resolution 


-0.0010 


-0.0003 


Event selection 






Multiplicity cut 


+0.0019 


+0.0026 
-0.0023 


Fiducial cut 


+0.0009 
-0.0028 


+0.0023 
-0.0024 


Monte Carlo statistics 


±0.0003 


±0.0007 


Total experimental error 


+0.0023 
-0.0032 


+0.0040 
-0.0042 
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Figure captions 



Fig. 1: The measured (a) EEC and (b) AEEC (points with error bars) compared 
with the JETSET (sohd histogram) and HERWIG (dashed histogram) Monte Carlo 
simulations. 

Fig. 2: EEC results: (a) SLD data fully corrected for detector and hadronization 
effects using JETSET (see text), (b) Detector correction factors from JETSET 
(solid line) and HERWIG (dashed line), (c) Hadronization correction factors from 
JETSET (solid hue with band) and HERWIG (dashed line). The width of the 
band in (c) represents the correction uncertainty calculated from JETSET (see 
text). Also shown as histograms in (a) are results of the fits of all four 0{ocl) 
calculations at / = 1 to the fully corrected EEC (see text). The fit range is 
indicated by vertical lines. 

Fig. 3: AEEC results: As Fig. 2, but (b) and (c) show the sizes of effective 
correction factors (see text). 

Fig. 4: (a) ag (M|) and (b) xlof Q^D fits to the EEC function 

of renormalization scale factor / (see text). The width of the band indicates the 
size of statistical errors, shown for the KN fit only. The vertical fines indicate the 
range used in the average (see text). 

Fig. 5: As Fig. 4, but for the AEEC. 

Fig. 6: Results of the fits of KN 0{al) and matched {modified R-matching) 
calculations to the fully corrected EEC. The dots represent the data points. The 
fit range is indicated by vertical lines. The solid histogram represents the fit of 
the matched calculation (/ = 1), the dashed histogram represents the fit of the 
0{a'g) calculation (/ = !), and the dotted histogram represents the fit of the 0{a^) 
calculation (/ = 0.0033). 

Fig. 7: (a) «s(A^|) and (b) Xdof ^^om fits to the EEC using 0{ag) and matched 
{modified R-matching) calculations (see text). 

Fig. 8: Results of the fits of the KN (!?(ag)+resummed calculation to the fully 
corrected EEC. The dots represent the data points. The fit ranges are indicated 
by the vertical lines. The solid histogram (dotted histogram) represents the fit at 
f — 1 using modified R-matching {R-matching). 

Fig. 9: (a) as{M^) and (b) Xdof ^om fits to the EEC using the resummed 
calculation combined with each of the four 0{ag) calculations using the modified 
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R-matching scheme. The width of the band indicates the size of statistical errors, 
shown for the KN fit only. The vertical lines indicate the range used in the average 
(see text). 

Fig. 10: As Fig. 9, but for the R-matching scheme. 
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